Neutralino dark matter and Higgs mediated lepton flavor violation in the minimal 

supersymmetric standard model 
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We re-examine the prospects for the detection of Higgs mediated lepton flavor violation at LHC, 
at a photon collider and in r decays such as r — > ^.rj, r — > /^7. We allow for the presence of a large, 
model independent, source of lepton flavor violation in the slepton mass matrix in the t ~ jj, sector 
by the mass insertion approximation and constrain the parameter space using the r LFV decays 
together with the _B-mesons physics observables, the anomalous magnetic moment of the muon and 
the dark matter relic density. We further impose the exclusion limit on spin-independent neutralino- 
nucleon scattering from CDMS and the CDF limits from direct search of the heavy neutral Higgs 
at the TEVATRON. We find rates probably too small to be observed at future experiments if 
models have to accommodate for the relic density measured by WMAP and explain the ((7 — 2)^ 
anomaly: better prospects are found if these two constraints are applied only as upper bounds. The 
spin-independent neutralino-nucleon cross section in the studied constrained parameter space is just 
below the present CDMS limit and the running XENON 100 experiment will cover the region of the 
parameter space where the lightest neutralino has large gaugino-higgsino mixing. 

PACS numbers: ll.30.Fs, ll.30.Pb, 12.60.Jv, 14.80.Ly, 14.80. Cp 



I. INTRODUCTION 

One of the appealing features of the minimal super- 
symmetric standard model (MSSM) with R-parity con- 
servation is the presence of a neutral, stable particle, the 
lightest neutralino, which presents all the characteristics 
to be a possible candidate for accounting for the cold dark 
matter in the Universe [l|. The amount of dark matter 
flh^, where fl is the dark matter density normalized to 
the critical density of the Universe and h is the reduced 
Hubble constant, recently has been precisely measured 
by the WMAP experiment 0. 

The Higgs sector of the MSSM Q, especially the heavy 
neutral Higgses A and H, play a prominent role in the 
physics of neutralino dark matter in two ways. In some 
region of the supersymmetric (SUSY) parameter space 
neutralinos yield the desired amount of relic density by 
annihilating into fermions through the s-channel reso- 
nant exchange of neutral Higgs bosons h, H, A, the so 
called funnel region where ttia — 2m^. As dark matter 
is expected to be distributed in a halo surrounding our 
galaxy, neutralinos can scatter off nuclei in terrestrial de- 
tectors: coherent scattering is mediated by scalar inter- 
actions through the s-channel exchange of squarks and 
t-channel exchange of the CP-even neutral Higgs bosons 
h and H. These effects become sizeable when squarks 
are heavy and tan /? is large in reason of the enhanced 
Higgs bosons coupling to down-type fermions, especially 
for the b quark which has the largest Yukawa coupling: 
moreover this couplings receive large radiative SUSY- 
QCD corrections at large tan /3 that can be relevant for 
their production in hadron-hadron collisions at TEVA- 
TRON and LHC In this scenario it is well known 
that S-mesons physics observables are very sensible to 
Higgs physics [ThiH and put strong constraints on the 



parameter space. The branching ratios for t lepton flavor 
violating decays are also enhanced near the experimental 
bounds [iMil, 0-126]. 

Once a source of LFV is present in the slepton mass 
matrix, for example the MSSM with the celebrated see 
saw mechanism for generation of small neutrino masses, 
two different mechanisms of LFV arise: gauge-mediated 
LFV effects through the exchange of gauginos and slep- 
tons [29l-[3l'| and Higgs- mediated LFV effects through ef- 
fective non-holomorphic Yukawa interactions for quarks 
and leptons [tI, [HI . LFV Yukawa couplings of the type 
L]^L-'j^H* are induced at loop level and become partic- 
ularly sizable at large tan /3. In this case the effective 
flavor-violating Yukawa interactions are described by the 
lagrangian: 



X (cos(/3 ■ 



cos^ /? 
a)h — sin(/3 



a)H-iA) + h.c. (1) 



where a is the mixing angle between the CP-even Higgs 
bosons h and H, A is the physical CP-odd boson, i,j 
are flavor indices that in the following are understood 
to be different {i ^ j). The coefhcients A*^ in Eq. ([T]) 
are induced at one loop level by the exchange of gaugi- 
nos and sleptons, provided a source of slepton mixing is 
present 

approximation are given by (21i] 



The expressions of in the mass insertion 
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where g and g' are the SU{2)l and U{1)y couphngs re- 
spectively, fj, is the the Higgs mixing parameter, mi^2 
the gaugino mass parameters and ni^^f^j^^ stands for the 
left-left (right-right) slepton mass matrix entry. /' is 
the derivative dl{x,y,z)/dz of the three point one-loop 
integral. The LFV mass insertions (5^^ — {m\y^ /ra\, 

^^RB. ~ (™lf,)*"'/"^?j' where {mj^ j^Y^ are the off-diagonal 
flavor changing entries of the slepton mass matrix, are 
free parameters which allow for a model independent 
study of LFV signals. 

The connection between gaugino-mediated LFV sig- 
nals and neutralino dark matter in the see-saw mecha- 
nism implemented in mSUGRA constrained MSSM has 
been recently studied in Refs. [l^, [2^: here it is shown 
that large neutrino Yukawa coupling affects the renor- 
malization group evolution equations of SUSY parame- 
ters from the grand unification (GUT) scale to the elec- 
troweak scale (in a SO(IO) GUT scenario) enhancing 
some LFV rates by orders of magnitude and changing 
also the neutralino relic density and direct and indirect 
detection rates. 

In this paper, on the other hand, we follow a differ- 
ent phenomenological approach: the study is done in the 
framework of MSSM with real parameters assigning the 
value of the parameters at the weak scale without any 
assumption on the mechanism of SUSY breaking or the 
high energy theory, nor on the origin of LFV entries in the 
slepton mass matrix and limitate our attention to Higgs 
mediated flavor violating effects. We study the inter- 
play between the assumptions of the lightest neutralino 
as dark matter candidate and Higgs mediated flavor vio- 
lation both to constrain the MSSM parameter space and 
to give prediction for the neutralino-nucleus scattering 
and LFV signals at colliders and in r decays. For related 
studies see |6l[27l[2i|. 

In Section|II]we discuss the scan of the parameter space 
in the real MSSM and the imposed constraints. Than we 
study their effects on the spin-independent neutralino- 
nucleon cross section and the arising correlations between 
supersymmetric parameters in Sec. IIIII In Section HVl we 
analyse LFV signals r — >■ /ury, r — > /i7, — >■ $ — ^ r/x -|- 
X at LHC and 77 — ^ r^bb at a future photon collider. 
Conclusions are given in Sec. [V] 

II. CONSTRAINED PARAMETER SPACE 

We introduce LFV in the model through the mass in- 
sertions Sj^j^ j^jzj = 0.5. This value ensures the largest 
rates in LFV processes and allow us to study the more 
optimistic scenarios; higher values contradict the mass 
insertion approximation as an expansion of propagators 
in these small parameters. Higgs mediated effects become 
eneteresting at large and tan/3 and low mA] further, 
if SUSY-QCD particles are heavy Higgs effects are dom- 



inant also for neutralino dark matter. We thus scan the 
following real MSSM parameter space: 

- 100 GeV <mA<l TeV; 

- 20< tan^ <60; 

- 500 GcV< ^ <5 TeV; 

The sign of /i is taken positive, as preferred by the 
SUSY explanation of the {g — 2)^ anomaly. 

- 150 GeV< mi, 7712 <1.5 TeV; 

We do not impose any relation but let them vary 
independently. To have large masses for gluinos we 
choose: 

1 TeV < TO3 < 5 TeV. 

- 1 TeV< m[/3 , m£)3, TTiQg <5 TeV; 

for the third generation of squarks: these are are 
varied freely without imposing any relation. For 
the first and the second generation the soft masses 
are set to be equal, = rnoi — ruQ. = rriq, 

where i = 1,2 and ruq is another free parameter 
which varies in the same range. 

- 300 GeV< mL^,mE^< 2.5 TeV; 

for sleptons of the third generation which are in- 
dependent parameters. For the first and the sec- 
ond generation the soft masses are set to be equal, 
rriL. — niEi — rrig, where i = 1,2 and is 
another free parameter which varies in the same 
range. Sleptons, diffrently from squarks, can be 
light in order to explain the (g — 2)^ anomaly. 

_2< "^'^3 <^ 2' 

— m[/3 ' m£)3 ' mE3 — ' 

for first and second generation the trilinear scalar 
couplings are set to zero. 

The outlined (16-dimensional) large parameter space 
is restricted imposing the following experimental limits: 

- Light Higgs and SUSY masses. 

LEP, TEVATRON bounds on sparticle masses and 
the LEP bound on light Higgs: > 114.4 

GeV [32]. 

- B-physics observables. 

For the i?-physics observables we use the the- 
oretical and experimental numbers of Table 1 
in Ref. [S^. Thus we require 0.995 < 

B«A-J(B^x 7) — 1-239, the upper bound on 
rare decays — > branching ratio is set 

to B{Bs ^ Ai+M') < 4.7 X 10-^ and 0.60 < 

^^sm < 1-24. Finally we require 0.85 < 

gA.s/Al(g^^^) ^^ 



There is at present a discrepancy between tlie SM value and the 
experimental value of the purely leptonic decays branching ratio 
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FIG. 1. Left: Scatter plot for the spin-independent neutralino-nucleon cross section versus the neutralino mass. The area above 
the sohd Une is excluded the CDMS final results; the area above the dotted line is excluded by the 2008 CDMS search. The 
dashed and dot-dashed lines give the sensitivity reach of two phases of the XENON experiment. The scanned parameter space 
and the imposed experimental constraints are described in Sec. |TI] The different graphical presentation of points corresponds 
to different steps in imposing the constraints on the {g — 2)^^/2 and on the relic density as explained in the legend of right 
panel. Right: Scatter plot in the {mA, tan /J) plane. The region delimited by the line is excluded by the CDF experiment. 



LFV T decays. 

Once Higgs mediated LFV effects are present in the 
model the non-observation of these rare decays puts 
strong constraints on the parameter space. The 
present experimental upper bounds are: B{t — )■ 
Hj) < 4.4 X 10-8 m, B{t ^ Air?) < 5 X IQ-^ H, 
B{t fififi) < 3.2 X 10-8 HI. 

Relic density. 

We use the conservative WMAP Str interval 0.09 < 
ilh"^ < 0.13 [1] on the relic density, both applying 
only the upper limit, (allowing for other sources of 
dark matter besides the neutralino) and the com- 
plete interval. See the legend of Figure [TJ 

Muon anomalous magnetic moment. 
The present discrepancy between a^^^ = {g — 
2')SM 12 g^jjj experimental measured value, 
Aa^ = a'^^'P - a^*^ lies in the interval (2 - 4) x 
10-9 [3i|. We always require al'^^^'' = {g - 
2^MSSM 12 < 4 X 10-9. We further show the mod- 
els which satisfy also the conservative limit lower 
bound aff^^^'' > 1 x 10"^^ ggg ^he legend of Fig- 
ure [TJ 



- Direct dark matter detection. 

The most stringent limit up to date in the 
neutralino mass range 100 — 1000 GeV for the 
neutralino-nucleus spin independent cross section 
comes from the CDMS experiment. The upper lim- 
its from the 2008 analysis [37[ and the recent final 
combined results [SSj are reported in Fig. [T] (left 
panel). 

- Non-standard Higgs search at TEVATRON. 
Recently CDF collaboration has published the most 
stringent exclusion limits in the {niA, tan/3) plane 
in the light of negative results in the search for 
heavy neutral Higgs bosons in the inclusive A, H 
production and the successive decay into t'^t~ 
pairs [39| . The excluded region is limited in the 
low niA, high tan/3 region and is depicted in Fig.[T] 
(right panel). 

For numerical computations we use the code Dark- 
SuSY [40] for accelerator bounds, the neutralino relic 
density and direct dark matter detection in our general 



B{B — > Tv) due to a recent analysis 13411 . Given the unclear 
situation both on the theoretical and experimental side we do 
not consider it here. 



MSSM. DarkSusy uses the code FeynHiggs [41[ for 
SUSY and Higgs mass spectrum and Higgs widths and 
branching ratios. For 6 — )■ 57 we used the routines in 
DarkSusy while for MSSM contribution to the muon 
anomalous magnetic moments {g — 2)^ those of Feyn- 
HiGGS which include also the leading and sub-leading 
two-loop contributions. For the others J3-physics ob- 
servales we used the formulas of Refs. [ifl, llil- We 
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generate 10^ random models, selecting the ones which 
evade the listed constraints. All of them are applied 
at the same time with the exception of the relic den- 
sity and (g — 2)^ anomaly for which we also relax the 
lower bounds: thus requiring only fl^h^ < 0.13 and 
a^issM < 4 X IQ-^ around 4 x lO'* survive, the light 
grey (turquoise) points in the Figures. Requiring also 
fi^/i^ > 0.09 around 7 x 10^ are left, the plus-shaped 
points, finally if aff^^^'^ > 1 x 10"^ only 52 remain, (the 
square points). 



III. NEUTRALINO DARK MATTER 

The spin-independent neutralino-nucleon cross section 
in the limit of heavy squarks and large tan /3 can be ap- 
proximated as d-d 



47rm^m^ 



tan^ 13 X Kf, 



(3) 



where iVn and N13 are the lightest neutralino unitary 
mixing matrix elements, rriN the nucleon mass (neglect- 
ing the mass difference between the neutron and the pro- 
ton) and Kf a factor which depends on nucleon form fac- 
tors. It scales like tan^ l3/m\ and it is able to constrain 
the low TO^-high tan/3 region even if to a lesser extent 
than flavor physics observables that scale like tan^ /3 / m\ . 

The right panel of Fig. [T] presents the allowed region 
in the (m^i, tan/3) plane: the region delimited by the 
line is excluded by CDF search in the channel pp — >■ 
A + X,A T+T~. The left panel of Fig. [1] shows the 
scatter plot for the spin-independent neutralino-nucleon 
cross section as a function of and the region excluded 
by CDMS [13, [11. We emphasize that CDF and CDMS 
limits are very mild constraints as can be seen in Fig. [1] 
The region excluded by CDF is practically excluded by 
the other constraints while the CDMS limit exclude only 
one plus-shaped point (not reported in Fig.[T]) leaving un- 
touched the regions preferred by WMAP and the — 
anomaly. Further, the final CDMS upper limits curve 
exclude around 300 light-gray (magneta) points between 
the solid and the dotted line in Fig. [T] and it is not still 
constraining the more interesting region. For clarity, in 
all the other plots only the final CDMS limits are ap- 
plied. Actually, it is more meaningful to compare the 
CDMS results with the plus-shaped and square points: 
in fact the limit on scattering with nuclei are extracted 
from rates which depend on the local density of dark 
matter in our galaxy halo which is assumed to be fur- 
nished by the weakly interacting massive particle, in this 
case the lightest neutralino. From this point of view the 
negative results of these experiments are natural in the 
present scenario and the two events found in the signal 
region by CDMS collaboration cannot be explained by 
our scenario. 

The XENONIOO experiment should reach the sen- 
sitivity corresponding to the dashed gray (red) line in the 



Figured] (left panel). Such sensitivity is able to cover the 
region with the highest cross section, > 300 GeV, 
where there is large higgsino component, as we will dis- 
cuss below. On the other hand the region preferred by 
{g — 2)^ anomaly cannot be covered. We also report the 
prospected sensitivity goal of the XENON experiment 
with 1 ton detector mass [42j, dot-dashed grey (red) line, 
which is 10~^^ — 10~^^ pb for neutralino mass in the range 
100 — 1000 GeV: practically all of the parameter space 
will be probed. 

As no relation has been imposed between the neu- 
tralino mass and tua and between gaugino mass param- 
eters, it is interesting to explore which correlations may 
emerge by the imposition of all the applied constraints. 

Fig. [2j left panel, presents the scatter plot of the ratio 
TO^/m^ versus mji/m^ where rriR is right-right mass pa- 
rameter for the stau. Points with the correct relic density 
abundance accumulate along the line mA/m^ ~ 2 where 
neutralino pair annihilation into fermions through reso- 
nant s-channel exchange of neutral Higgs bosons A, H 
is the dominant mechanism in large portion of the pa- 
rameter space. Stau coannihilation is at work for models 
where mn ^ m-^ and coannihilation with the second neu- 
tralino and the lightest chargino are important for larger 
values of the ratios. 

In the right panel of Fig. [2] the gaugino fraction 
|A^iiP + |A^i2p is plotted against the neutralino mass. For 
neutralino masses below 400 GeV, the preferred region 
by the (17 — 2)^, is pure gaugino while for masses greater 
than 400 GeV higgsino component is present. This effect 
can be seen in Fig. [T] in the spin-independent neutralino 
cross section which depends on Ni^ (|A^i3p + \Ni 



1 



|A^i2p): the models with the highest cross 



section are the ones with > 400 GeV where the xx^ 
coupling is enhanced in reason of a larger higgsino com- 
ponent. 

The left and right panels of Fig. [3] present the scatter 
plot in the (/i/mi, 1112 /mi) and (/i/m2, rfi2lrfi\) planes 
respectively. We see that models with the correct relic 
abundance have TO2 > mix /i > m\. The (5 — 2)^ prefer 
models with /i > m^. The models with strong gaugino- 
higgsino mixing, mi ~ /i, ^ < m2, |A^iiP + |A^i2p < 0.9 
can be probed by XENONIOO. We further note that 
most of the point in WMAP and {g — 2)^ ranges are 
charactherized by having high degeneracy in the gaug- 
ino masses m\ ~ TO2. Such conditions give a "well- 
tempered bino/wino" neutralino which can be realized 
in model with split supersymmetry as shown for exam- 
ple in Ref. [43|] 



IV. PROSPECTS FOR LEV SIGNALS 

The T LFV t decay which is more sensible to Higgs 
mediated effects is t ^ fiij [13,, J^, ,24-,26i] and the 
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FIG. 2. Left: Scatter plot of the ratio m a /m^^. versus the ratio m_R/m^. Right: Scatter plot of the gaugino fraction |Afiip + |Afi2p 
versus the neutralino mass. The scanned parameter space and the imposed experimental constraints are described in Sec. |TT1 
The same legend of Fig. [T] applies. 




FIG. 3. Left: Scatter plot of the ratio fi/rrii versus m2/m\. Right: Scatter plot of the ratio iJi./m2 versus m2/mi. The scanned 
parameter space and the imposed experimental constraints are described in Sec. |lTl The same legend of Fig. [T] applies. 



branching ratio reads [la, |2C 



97r^ 



tan^ (3. 



(4) 

Here ~ 110 MeV, B(t ^iw) ^ ^^d F,, 

a factor which depends on the hadronisation of quarks 
bilinears matrix elements: in the treatment of Rcf. [IG'I 



it is such that F,^ 
through the factor 



20 Lepton flavor violation enters 



(5) 



The approach using chiral perturbation theory llTI . |26(| gives re- 
sults different at most by a factor two. These uncertainties will 
not change our conclusions. 
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FIG. 4. Left: Scatter plot of the LFV violating vertex in Eq. (O as a function of tua- Right: Scatter plot of the branching 
ratio $ — > r/i versus tua ($ = A, H). The scanned parameter space and the imposed experimental constraints are described 
in Sec. [Ill The same legend of Fig. [T] applies. 



We also consider the radiative decay r — >■ /i7 which re- 
ceives also important contributions by gaugino-slepton 
loop diagrams: the factors |A|^p and |A|j^p enter sepa- 
rately in the branching ratio and not through the com- 
bination in Eq. (EL For the computation we used the 
formulas in Ref. [21| including both gaugino mediated 
and Higgs mediated effects. 

The CP-odd Higgs boson decay width and branching 
ratio are [Hi 



i tan^ ^ A^ r(A ^ r+T 
tan"^ (3 A^B {A ^T+T 



(6) 



where A^ is defined in Eq. ([5]) and we used the fact that 
r{A r+/i-) = T{A r-/i+). For the CP-even Higgs 
boson H, the right hand sides of Eq. ([6]), are multiplied 
by a factor sin(/3 — a)/(cosa)^ which is order one in our 
scenario where tua — ttih, thus the previous formulas 
hold for both bosons. 

In Fig. IIVI left panel, we present the scatter plot for 
the effective vertex given by Eq. ([SJ , while the right panel 
shows the scatter plot of the branching ratio given by 
Eq. as a function of to^. We see that they reach 
10~^ and 3 x 10~^ respectively for models preferred by 
WMAP measurements, two orders of magnitude lower 
than what found without imposing it on the parameter 
space [i^, m, 11^ . The branching ratios of LFV decays 
T — ^ /i?7 and t — >■ /i7 versus are given in Fig. [5j For 
models with the correct neutralino relic density abun- 
dance and preferred by the {g — 2)^ anomaly, both 
are generally under 10"^'^ while relaxing the constraints 
lower bound they can reach the 10~^ level. We remind 



that at a Super-B factory the present limits 0(10^*) can 
be lowered to O(10~^ ~ 10~^°) for the ixrj hnal state 
because the branching ratio scales linearly with the lu- 
minosity due to practically negligible background. In the 
/i7 case for the presence of large background the branch- 
ing ratio scale as the square root of the luminosity and 
the sensitivity reach is one order of magnitude lower j4jj . 

We further revisit the prospects for detection of Higgs 
LFV signals in pp collisions at LHC [l^ HI] and in 77 
collisions at the photon collider option of the future In- 
ternational Linear Collider [46l |. 

At high tan (3 the dominant production mechanisms for 
A, H at LHC is hh fusion due to the mb tan (3 enhanced 
66$ couplings. We calculate the cross section with Feyn- 
HiGGS which uses the approximation 



MSSM 



{bb ^ $) = cr^*^ (66 ^ $ 



66) 



MSSM 



r($ ^ 66) 



SM 



(7) 

where a^^'{bb — ^ $) is the total SM cross section for 
production of Higgs boson with mass m$ via 66 fusion: 
to obtain the value in the MSSM it is rescaled with the 
ratio of the decay width of the inver se p rocess in the 
MSSM over the SM decay width mHliil. We calculate 
for each random model the product the cr{pp — $ -|- 
X) X S($ Tfi). As masses and couplings of A and 
H are practically identical as discussed above, we have 
a{pp A + X)+ a{pp H + X)~ 2a{pp ^ A + X). 

The scatter plot a{pp — > $ + X) x B{^ Tfi) is shown 
in Fig. ini left panel. We see that with the nominal in- 
tegrated luminosity of 100 fb~^ per year models which 
satisfy both the relic density abundance and Aa^ can 
give up to 10 events per year (squared points), up to 40 
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FIG. 5. Left: Scatter plot of the LFV violating vertex in Eq. (O as a function of tua- Right: Scatter plot of the branching 
ratio $ — > r/i versus tua ($ = A, H). The scanned parameter space and the imposed experimental constraints are described 
in Sec. [Ill The same legend of Fig. [T] applies. 



if we relax the condition on the lower limit of Aa^ (plus- 
shaped points) and up to 200-300 relaxing both the lower 
limits (magenta (grey) points). 

In 77 collisions the main production mechanism for 
^ ^ A, H is TT fusion [l^l while the bb is suppressed 
by a factor 3(l/3)^(TOh/TOT-)^ ~ 0.1 which cannot be 
compensated by corrections to the b Yukawa coupling. 
In Ref. [4^ we studied in detail the fir fusion process 
77 — iirbb where the Higgs boson is produced in the 
s-channel via a virtual /ir pair and can be detected from 
its decay mode A — >• bb. 

A good analytical approximation for the cross section 
is obtained using the equivalent particle approximation 
wherein the colliding real photons split respectively into 
T and /i pairs with the subsequent /zr fusion into the 
Higgs boson. The splitting functions of the photon at 
leading order read [47l] : 



thus the cross section is given by: 

47r2 r(^ ^ Tfi)B{A bb) 



(8) 



ct(77 fiTbb; s^^) 

- In l/t 



(9) 



with t = iTiA/'iE-f, ry = In ^Jx^/xr, x is the fraction of the 
energy of the photon carried by the virtual lepton. In (46j | 
we have shown that the effect of photons spectra can be 
neglected, we thus consider monochromatic photons with 



= 600 GeV, and photon-photon luminosity 500 
based on the parameters of TESLA(800) 
The scatter plot of the signal cross section versus ruA is 
shown in Fig. |51 right panel. Here the models which sat- 
isfy both the relic density abundance and Aa^ (squared 
points) have maximal cross section 10""^ fb, which is too 
small. Relaxing the lower limits cross section values up 
to 2 X 10"^ fb are possible, giving 10 events/year. 



V. SUMMARY AND CONCLUSIONS 

In the framework of the MSSM with heavy SUSY-QCD 
particles and large tan (3 we have studied lepton flavor vi- 
olation mediated by the heavy neutral Higgs $ = A, if in 
T~ ^ sector both in t lepton decays, r — >■ /i?] and r — >■ /X7, 
and at high energy collider through the production and 
decay at LHC, pp ^ A- X, ^ ^ and the ^ — r fu- 
sion at a photon collider, 77 — > r^bb. The approach to 
LFV has been model independent by the use of the mass 
insertion approximation. We used large mass insertions 



"ll 



^RB. ~ estimate the number of events in 



the most favourable scenario that can be obtained in fu- 
ture experiments. With such a large source of LFV the 
branching ratios of rare decays can exceed the present 
experimental upper bounds from BABAR and BELLE 
which therefore provide constraints on the MSSM pa- 
rameter space in presence of LFV. Other constraints that 
have been imposed are limits from direct search of spar- 
ticles and of the light Higgs, i?-physics observables, the 
{g — 2)^ anomaly, and recent limits from TEVATRON 
search of non standard Higgs bosons in the tt channel. In 
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FIG. 6. Left: Scatter plot of the inclusive production cross section pp ^ $ + X times the branching ratio of <1> — >■ r/i at 
LHC versus ttia = A, H). Right: Scatter plot for the cross section of the process 77 — > r^bb in photon-photon collision at 
ys-y^ = 600 GeV. The scanned parameter space and the imposed experimental constraints are described in Sec. [Ill The same 
legend of Fig. [1] applies. 



the R-parity conserving MSSM the heavy neutral Higgs 
play an important role both in neutralino annihilation 
cross sections to satisfy the relic density of dark matter 
measured by WMAP and in the spin-independent neu- 
tralino nucleon cross section in direct dark matter search 
experiment. We thus have imposed on the MSSM pa- 
rameter space the present CDMS exclusion limit and the 
WMAP limits on ilK^ . 

We have found that in models with 0.09 < rj^/i^ < 0.13 
and 1 < aff^^'^'^ x 10^ < 4: the branching ratios of r — >■ 
r — > //7 are both under the 0(10""'^°) thus probably 
undetectable even at super-B factory; at LHC the cross 
section for pp — )• ^+X, $ — t/z can reach O(10~^ — 10~^) 
fb in the range tua = 400 — 1000 GeV giving up to 10 
events with 100 fb~^; the cross section of 77 — >■ r/i66 
reach 0(10"^) fb, thus too small even for the large value 
of the expected luminosity of 500 fb^^. Prospects are 
somewhat more encouraging if we relax the lower limits, 
imposing only fl^h^ < 0.13 and aff'^'^'^'^ x 10^ < 4. In this 
case branching ratios of LFV r decays can reach 0(10^), 
the cross sections at LHC about 2 fb for low itia masses 
and around 2xlO~^fbin77 collisions. 

We derive two main indications from this analysis. On 



one hand, even with large sources of lepton flavor viola- 
tion in the slepton mass matrix, the process where Higgs 
mediated t — /i transitions should manifest could be be- 
yond the sensitivity reach of future experiments. On the 
other hand, to observe such effects, in any case, the full 
luminosity of the machine is needed. 

We emphasize that our not optimistic conclusions are 
specific to Higgs mediated effects. As shown in [23|, in 
typical SUSY parameter space where gaugino-mediated 
effects are dominant over Higgs mediated ones and in 
the context of SUSY see-saw mechanism, LFV rates are 
detectable by future experiments. 

We have also studied the spin-independent neutralino 
nucleus cross section: we have shown that in models that 
satisfy 0.09 < n^h^ < 0.13 and 1 < aff'^'^^ x 10^ < 
4, the cross section lies just below the sensitivity of 
XENONIOO which should report results soon. The full 
XENON 1 ton is needed to cover all the parameter space. 
However, if the lower limit on (5 — 2)^ is not considered 
XENONIOO is sensitive to the neutralino mass range 300- 
1000 GeV in models where the higgsino component is 
large. 
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